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Abstract We show by explicit computation that the recently discovered duality invariance of D=4 
linearized gravity fails, already at first self-interacting, cubic, approximation of GR. In contrast, 
the cubic Yang-Mills correction to Maxwell does admit a simple deformed duality. 
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! 1 Introduction 

O 
IT) 

"Duality" has become a touchstone in relating seemingly different regimes and models in field and 
string theory. Its humble origins lie in the ancient observation that in (and only in) D=4, on-shell 
configurations with mutually rotated electric-magnetic fields also obey the source-free Maxwell 
equations. While this observation is essentially correct, it is sometimes misunderstood. Duality 
can only be properly formulated in terms of the unconstrained dynamical variables, rather than 
as a formal E <-> B rotation [1]. This is already apparent from the facts that B is identically 
diver genceless, and so can only be related to E upon implementing the latter's Gauss constraint, 
and from the Maxwell Lagrangian's second order, "hyperbolic", non-dual invariant, form Lm = 



X 



Ue 2 -b 2 ). 



Surprisingly, duality invariance was recently [2] extended to free massless spin 2, then to 
free gauge fields of any spin and statistics [3]. Can this abelian duality invariance be promoted 
to encompass the two physical self-interacting generalizations, vector (YM) and tensor (GR)? [No 
consistent nonabelian higher spin models are known.] The YM conjecture was already considered, 
and settled in the negative in [1], although we will see that, actually, its first-cubic-extension does 
allow a natural, consistent deformation of duality, at least in Coulomb gauge. Our other, and 
principal, objective is to settle the question (raised in [2]) of extending duality to GR. We will show 
(somewhat laboriously) that this cannot be achieved, at least perturbatively: (deformed) duality 
already fails at its first, interesting, cubic, level. 

We will begin with a brief summary of what linear duality is (as well as of what it is not) for 
free gauge fields: the two degrees of freedom of all D=4 free gauge fields-their ±s helicities-can be 
rotated into one another, by a canonical transformation mixing their two pairs of unconstrained 
dynamical variables, while keeping the Hamiltonian form-invariant. This is not to be confused with 
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other transformations, such as the "harmonic oscillator" (p — ► q, q — ► —p) rotations within a given 
mode. We will then revisit the YM system and establish its cubic order duality invariance, before 
turning to the GR case. An Appendix provides details of the latter. 



2 Free Vector Field Duality 

Free gauge field duality is a canonical transformation, linking coordinates and momenta of different 
excitations, that leaves the Hamiltonian form-invariant. The simplest case is Maxwell's, whose 
first-order action is 

W£ r ,A r ]=/d 4 *[(-i^ (2.1) 
J J LA=1 

upon implementing the Gauss constraint, V ■ E = 0. Here, transverse vectors are labelled with a 
"T" , from the usual transverse-longitudinal decomposition, 

V = V T + V L , V -V T = = VxV L , J d 3 xV T -W L = . (2.2) 

The indicated orthogonality between any pair of T and L vectors implies that only the manifestly 
gauge-invariant A T component survives in the action (2.1). The duality rotation's infinitesimal 
form is (we exhibit the - redundant - T-index on B for emphasis): 

5E T = B T = V x A T , 5A T = V~ 2 V x E T =► SB T = -E T . (2.3) 

This E T <-> B T rotation clearly leaves / {E\ + B\) invariant. The middle equation demonstrates 
that the desired B-rotation is indeed implementable at the level of the canonical coordinates A T . 
That (2.3) is also canonical, i.e., that the symplectic form S = J pq term is invariant, follows 
from the fact that the curl and Laplacian are hermitian operators, O, by virtue of which any 
/ X@X = — J 'XPX = — I X@X = upon (double) parts integration. The loss of manifest Lorentz 
invariance and of space locality inherent in this procedure is entirely harmless and indeed necessary 
even to formulate, let alone establish the transformations as canonical ones. [Note that it is only 
in three space dimensions that the vector E can even be matched with the magnetic tensor F^j (by 
dualization in the e* jfc sense).] 

The above rotation is quite different from "harmonic oscillator" duality, valid in any dimen- 
sion 

I[m, k]= J\pq- \{p 2 jm + kq 2 )] , 

Sp = aq , 5q = —(3p , a/(3 = km . (2.4) 

that relates a single excitation's variables, and reflects the equivalence of different parameter 
regimes 1 through the dependence of the solutions on (k/m). There is a similar "duality" of Maxwell 
theory: rewriting (2.1) as 

I M = J {(-E).A-\[E 2 + A(-V 2 )A)} , (2.5) 

immediately implies invariance, within each helicity, under 

5E = V^V 2 ~A, 5A = -{-V 2 )- 1/2 E . 

x Or relating different models, here I[m,k] <-> 7[fe _1 ,m _1 ], with different parameters [4]. 
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3 Cubic Yang-Mills Duality 



As an instructive (and transparent) contrast to GR, we study first the nonlinear extension of a 
Maxwell multiplet by adding the cubic terms in YM, and show that, surprisingly, they permit-a 
deformed version of-abelian duality. [This in no way contradicts the demonstration in [1] that full 
YM precludes duality.] For simplicity, we work with SU2, whose structure constants e abc permit an 
obvious 3D internal vector notation. The first-order covariant YM action is (setting g = 1), 

Iym[F,A] = -\J [F^ -(d^A u -d v A^A^A u ) - \F^-F^ V ) . (3.1) 

The (3+1) versions of (3.1) and of the Gauss constraint become 

Iym[E, A] = \ J[(-E)-A-\{E 2 + B 2 )], B = VxA + AxA, DE = (V+Ax) -E = 0. (3.2) 

Adopting Coulomb gauge, A L = simplifies the process: 

Iym = \\ {(-E t )-A t -\[E 2 t + E 2 l + B 2 )) , (3.3) 



(V + A T x) ■ E L = -A T xE T ■ (3.4) 

The constraint 2 fixes El in terms of the dynamical pairs: El = + O(A T E T ); since E 2 L is quartic, 
we may drop it from (3.3) to cubic order; omitting "T", this leaves 

Iym[E,A] = I M [E,A]+I YM [E,A] * \ J {(-£)-A-i(£ 2 +(VxA) 2 } + J VxA-AxA . (3.5) 

Only the final term in H differentiates the action from that of a triplet of photons. The original 
linear duality transformation (2.3) clearly alters (only) this term, 

S L J(VxA-AxA) =3 J E-AxA^O . (3.6) 

To cancel this cubic term, we must deform the abelian transformation by adding a quadratic SqE, 
that will generate a cubic variation from the / E 2 in (3.5); the obvious choice is 

5 Q E= (AxAf . (3.7) 

[We have projected the "T" part of (Ax A) since E is transverse; this is just a formality here, and 
throughout, since orthogonality would automatically perform the projection in / E -SqE.] Having 
succeeded in keeping Hym invariant, we must still check that the symplectic variation, generated 
by (3.7), namely 

Sq j EA = j AxA-A (3.8) 

vanishes. Indeed, time integration by parts shows it to equal minus twice itself. 

To summarize, we have succeeded in keeping duality invariance of YM to lowest nonlinear 
order, in Al = gauge at least. This was accomplished by setting 5E = B = 'V x A + Ax A, 



"Formally, inversion of (3.4) for El[At, Et] would solve classical YM entirely! 
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certainly the most obvious guess; the corresponding part of the generator is just the YM Chern- 
Simons form: 

G E = j d 3 x[A-VxA + lA-AxA] , 

5E = [G E ,E] = VxA+AxA = B YM ■ (3.9) 

We have not attempted to extend this process to quartic order, where terms ~ (A x A) 2 from B 2 
and ~ (A ■ E?) 2 from E 2 L appear, nor do we know whether the cubic order success has some deeper 
physical origin. 



4 GR Duality Fails 

Having illustrated how duality deformations can succeed at cubic level for YM vectors, we now turn 
to the-considerably more complicated-tensor case. We will begin with a description of Pauli-Fierz 
free spin 2 theory, transcribed into a notation manifesting its duality invariance. We will then derive 
the cubic correction to GR, and subject it to the abelian transformations. Finally, we will show 
that its lack of invariance under the latter cannot be compensated by adding further, quadratic, 
deformations of the dynamical variables: there is no (perturbative) duality invariance in GR. 

Let us first express the familiar free spin 2 gauge system in "dual-ready" form, using the first 
order formulation of full GR [5] , thereby also obtaining the cubic correction in a unified way. The 
Einstein action is 

I E [7r,g] = I d 4 x [^ 9ij - 

N = (-g 00 )- 1 / 2 , Ni = goi 

R° = - \fi~g 3 R+ - 5*i*i)( 3 g)- 1/2 , /." 2 . (4.1) 

All explicit (roman) indices refer to the intrinsic 3-space, of which 3 i? is the scalar curvature and Dj 
the covariant 3-derivative; it 13 , essentially the second fundamental form (density), is an independent 
variable on a par with g^. Our convention is Rij ~ +8^^. The expansion about flat space is 3 

hij = gij -Sij , N = 1 + n ; (4.2) 

instead, Ni and it 13 vanish in flat space and are of first order. Abelian gauge invariance and use 
of the four constraints, R^ = simplifies (4.1) dramatically. We recall the familiar transverse- 
traceless orthogonal decomposition of a flat space (where index position is immaterial) spatial 
tensor = Tji, 

Tij = T^ + ^-V^d^TT + ^ + djTi) 
d t T? T = = If T , djiSij - V" 2 4) = ee (S^ - V 2 d 2 )d 3 . (4.3) 

Of the six components of , two are "TT" one is "T" and three represent the vector T • Orthogo- 
nality between the various components under integration is manifest; one example is / d s xT^ T diWj = 

3 We spell out the fact that gij is exactly Sij + hij to all orders, so that neither hij nor 7r lJ are subject to further 
expansion, unlike say g %i ~ Sij — h lj + 0(h 2 ) or y^g ~ 1 + \ ha + 0(h 2 ). 
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0. Once decomposed in this fashion, the twelve (n lJ ,hij) components are easily classified: At lin- 
earized level, the theory's abelian invariance under 5h^ u = [d^ v + d u £a), simply means that the 
four gauge variables (ir , hi) do not appear in the quadratic action. Likewise, the four linearized 
constraints, 

-,R -► ^L L {h) = (dfjhij - V 2 hu) = -V 2 h T = , K -23^ = (4.4) 

remove the respective components (h T ,7r t ), as is also obvious for ir 1 upon using (4.3). The four 
remaining (tt ij , hij) then simply reduce to their two pairs of "TT" degrees of freedom. We designate 
them by 4 (p*- 7 , hij) henceforth dropping the "TT" notation. The symplectic term is just that of the 
two excitations, 

S = J <l\r - J d 4 xp ij qij , (4.5) 

while the Hamiltonian is simply the quadratic part of R°, subject to the linear constraints. A short 
calculation (essentially integrations), finally reduces (4.1) to the usual harmonic oscillator form, 

ii \P, q]=jd 4 x^ PAQ A ~W + ( Vg) 2 ] | , (4.6) 

with traces on implicit indices understood. For our purposes, the notion of duality in this symmetric 
tensor world emerges from the existence of a generalized curl [3], distributed on the tensor's indices: 

(OT)^ ee ±(e am d e T mj + e^ m d,T im ) = (OT) Jt . (4.7a) 

Its nonlocal extension 

(07% ee (V- 2 07% (4.7b) 

will also streamline notation; both 0-operations are hemitian. Thus, acting on the curl defines 
a "magnetic" field, 

Bij = {Oq)ij = Bji , (4.8) 

in terms of which (4.6) becomes 

1%, q]=J d 4 x {pq - i(p 2 + B 2 )} . (4.9) 

This form makes it irresistible to perform the rotation 5 

Spij = B,^ , Sqij = (Op)ij 5Bij = -pij , (4.10) 

the final step uses OO = —V 2 , just as V x Vx = —V 2 on transverse vectors. This transformation 
obviously maintains the Hamiltonian of (4.6), while the g-sector of (4.10) shows that 5B is indeed 
implementable as a transformation of the underlying coordinate q. Invariance of the symplectic 
term is guaranteed (as for spin 1) by hermiticity of the O and O operators. 

Establishing duality transformation and invariance of free spin 2 under them was the easy 
part; now we must return to the full action (4.1) to obtain its cubic, I%\p, q], correction, subject it 

4 More precisely, the original linear Hamiltonian is + -j(Vftij) 2 ], requiring the rescaling \f2-K^ = p IJ , 

l/\/2 hij — qij to bring it into standard free-field form (4.6). 

5 Higher free spin fields [3] follow the same pattern, in terms of suitably generalized curl operations to accommodate 
any number of spatial indices. 
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to the linear rotations (4.10), and try to cancel away their (nonvanishing) effect by quadratically 
deforming them. 

The cubic correction to H is formally easy to find. It is simply the cubic part, R®, subject 
to the lower level constraints. The procedure is greatly simplified, as it was for YM, by choosing 
"Coulomb" gauge, here hi = = ir . [Recall that only the linearized action is invariant under 
the linear part of the non-abelian 8g^ v = D^ v + The upshot of this process is the cubic 

Hamiltonian 

H c (p,q)=R° c (p,q) = J d 3 x 

the last term in (4.11) vanishes due to p im Se m = 0. Furthermore, the symplectic term has the great 
virtue of remaining quadratic in our gauge, while the lower constraints also imply that the same 
(p, q) are dynamical. 

Whereas the cubic order curvature scalar density is a bit elaborate, we are in fact only 
interested in its (linear) variation; 

h J ( 3 RV9)c = ~ J (VgG ij ) Q b L q l3 , (4.12) 

by the Palatini identity. The variation of the cubic kinetic terms is just 5l J 2trppq, so we have 

S L H c (p, q)= J d 3 r GpLpV + 2 J (pp^LQij + 2 J {6 L p lj (PQh + (php)^} ■ (4.13) 



3 R^g) c + 2p ij p em (q ii S jm - \qijStm) 



(4.11) 



Schematically, the variation (4.13) is of the form 5lI ~ q 2 p+p 3 . There is an a priori daunting 
array of possible compensating quadratic deformations of each (p,q): o~Q(p,q) ~ (p 2 + pq + q 2 ), 
with various possible derivative and index structures in each. Further, all are subject to keeping 
the relative dimensions of {p ~ q, q) as well as enforcing 5 J pq = 0. Owing to the dimensional 
asymmetry between 5qp and 5qq, the latter will involve the non-local operator 0, as compared to 
Sqp ~ O. This will lead to the presence of both local and nonlocal variations that must separately 
cancel. 

[Parenthetically, we lay to rest the otherwise attractive idea that, since duality invariance 
requires the other Poincare generators of the cubic model to be invariant too, why not test the 
momentum generator instead, say? Unfortunately, Dirac's dictum that the Hamiltonian form 
of dynamics always has a simple momentum holds here as well. The transformation generator 
G = Pi5x % , with Pi = —2 j dj-rrj and the constraint Dj-ir 1 ^ = says that to all orders, 

djirj = Yj\ n{ = [ji, k]^ k = \pV t q . (4.14) 

Hence P = J pVq to all orders [5] and its invariance provides no independent test beyond that of 
/ PlftQ- T ne same holds for the rotation generator, of course.] 

In the following, we will just outline the flow of possible deformations and their consequences. 
The more unpleasant details are relegated to the Appendix. Returning to (4.13), we have 

d L H c (p, q) = J d 3 r [G% (Op)ij + 2tr ppOp + 2tr {qpOq + qOqp}] 

~ p 3 +q 2 p, (4.15) 
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Obviously, the simplest term to cancel is the pure p 3 , which can (only) be accomplished by S^qP ~ 
p 2 . More precisely, acting on the quadratic Hamiltonian, this deformation leads to 6 

S^pn ~ {p6p) tj <5g } Jp 2 ~ tr(pp6p) , (4.16) 

exactly of the form of the term to be cancelled. Now, however, we have to cancel the unwanted 
effect of (4.16) on / pq, 

Jpq-trJ {p6p) tj q i3 . (4.17) 

The unique possible cure for this is 

S^q-iqOp+pOq)*,; (4.18) 

and indeed it works, cancelling (4.17). While we have now gotten rid of the p 3 part of (4.15), we 
must not forget one further effect, that of (4.18) on Hq, 

S Q1 J %'(-V 2 )^ ~trj q{-V 2 )[qdp + pdq\ ~ q 2 p , (4.19) 

a generically nonlocal term. The possible remedies to the overall / q 2 p residue of (4.15) plus (4.19) 
consist of Sq^p ~ q 2 and 5q\ ~ qp. The latter choice is clearly not desirable, since we have just 
exploited it in (4.18), in an essentially unique way. We are therefore stuck with the last hope, 

6$>p ~ qOq , 8% \Jp 2 =tr J p(qOq) . (4.20) 

Now we exploit the local/nonlocal division of the variations. In (4.15) there are both kinds; the 
former type is easy, since (by dimensions) it has no explicit derivative beyond the curl O, that is, 

SiocalHc ~ 2 J {p(Oq)q + qOqp) . (4.21) 

This is clearly cancelled by (and only by) 

S^p ~ A{qOq) , 5\Jp 2 = ttr J p{Oq)q . (4.22) 

The symplectic contribution of (4.22) "miraculously" vanishes; its (local) form is 

€lljp'i-^l( q q) ij (Oq) ij . (4.23) 

The positions of the dot and the curl are immaterial. What matters is that O is hermitian, while 
d/dt changes sign on integration by parts, which suffices to show this term is proportional to minus 
itself and vanishes. The (decisive) nonlocal terms require the Appendix! 

6 Here, and in all quadratic transformations, a projection to "TT" space is understood, to keep the character of 
(p, q). However, this is really immaterial, since they always safely multiply a linear TT variable. 
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5 Summary 



We have considered whether any deformations of the duality transformations of linearized vector 
and tensor gauge theories might rescue (an extended version of) duality in their nonlinear regimes. 7 

For YM, it was surprisingly easy, at least in Coulomb gauge, to do so at leading, cubic, "post- 
Maxwell" order, simply by letting the electric field variable rotate into the full YM magnetic field, 
while keeping the vector potential's rotation unaltered. We have not explicitly analyzed quartic 
order and beyond in the theory's infinite series expansion in (p, q) , since it is known that duality 
fails for full YM. 

Extension of free spin 2 duality to cubic order proved considerably more complicated to 
decide, but we were able to show that no deformation compensates for the loss of abelian duality 
invariance, even in Coulomb gauge. We conclude that, perhaps disappointingly, rotation among 
helicities ceases to be an invariance beyond the free spin 2 level. 

This work was supported in part by NSF grant PHY04-01667. 



A Appendix 



In section 4, we sketched problems of extending linear spin 2 duality to the cubic term of GR. 
This Appendix provides details of the obstacles and shows why they cannot be overcome. We first 
collect the key quantities, starting with the gravity Hamiltonian at quadratic 



He 



d 3 x 



(Vq) 2 +p 2 



and cubic, 



H c (p,q) = J 



i Rc(q) + 2p t n jf m q it 



(A.l) 



(A.2) 



levels, with all relevant constraints satisfied: only "TT" variables are involved. [The corrections 
coming from solving the nonlinear constraints affect only quartic terms.] The abelian duality 
transformation (4.10), which preserves Hq, is 



6^q ab = (6p) ab = P ab 



5^p ab = (Oq) ab = Q ab . 



(A.3) 



We now consider the effect of this rotation on H c , starting with the variation of the simple cubic 
vertex ppq of (A.2) 

-25 j ppq = -2 J 5 (1 kabP ac p b c + 2q ab 5^p ac p b c = A + B. (A.4) 
The A term's contribution in the above variation is proportional to p s . Specifically we have 



-2(5^q ab ) ^ b c = P ab (pp) ab . 



(A.5) 



7 The vector and tensor models are of course very different, not only in index complication, but in the fact that 
cubic and higher YM terms are of lower derivative order, whereas in GR, all order have two derivatives. 
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The unique and obvious way to cancel the above variation introduces a contribution in 5p quadratic 
in momentum, 

S^pbc = -PabPbc ■ (A.6) 

The ensuing variation of the p 2 term in Hq then exactly cancels (A. 5). This deformation of 5p 
obviously affects the symplectic form S, which in turn must be compensated. Indeed, 

5S = J d 4 x 5^p ab q ab = - J d 4 xP abPbc q ac . (A.7) 

Cancelling this variation requires, and completely fixes, a quadratic term ~ pq in Sq ab - [Any 
alteration would automatically produce a change in the p 2 term 5^p ab , which, however, has been 
just determined.] To evaluate the required qp term in q ab , we first integrate (A.7) by parts with 
respect to time, 

5S = J d A x tr[qPp + (pq)P] . (A.8) 
Thus, the desired variation must be 

5$q ac = [q cb P ba + (a&)] (A.9) 
Next we consider the B term in (A. 4), namely 

-4 J d 3 x q ab (5 (1 V C ) p b c = -Aj d 3 x tr(qQp) . (A. 10) 

Variations of the same pq 2 form will also contribute from 3 R c {q) in (A. 2), 

-5 J d 3 x^R =-J d 3 x^G ab {q)S (1) q ab = 2jd 3 x R$ P ab , (A.ll) 

as well as from the q 2 term in Hq 

-\jd 3 x q a bV 2 5 (2 }q ab = ~hjd 3 x tr[qV 2 qP + qQp] . (A. 12) 

The variation in (A.ll) can be re-arranged by explicitly displaying the Ricci tensor 8 to 0(q 2 ): 

^R ab = (I q ac V 2 q c b + \ q bc V 2 qi) - \ d r (q rn (d a q nb + d b q na - d n q ab )) + \ d a d b (q rn q rn ) 

- \ d a q mr d b q mr - \ (d r qr - d m ql) (d m q ar - d r q am ) . (A.13) 

The 8 a 8 b term can be dropped because it is longitudinal and vanishes when integrating with a TT 
object. The variation (A.ll) can be now rewritten as follows 



q ac ^ 2 qt + qbc^ 2 q c a - di(q ln {d a q nb + d b q na - d n q ab )) 



2 J d 3 xR ab (q)P ab = J d 
-\ d a q m id b q ml - (d m q a id l q bm - d m q a id m q l b ^ 



Pab ■ (A.14) 



8 We are allowed to drop (g ab R) Q = 8 ab R Q — h a bR.L; the first of these vanishes when contracted with a TT tensor, 
the second because Rl = is the linear constraint. 
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The combination of the variation (A. 14) and (A. 12) simplifies a bit to 

J d 3 x \{-di(q ln {d a q nb + d b q na - d n q ah )) - \d a q mi d h q ml - d m q a id l g bm | P ( 

+ \{qa C qb)Pab - \qTQabPa; 

Finally, inclusion of (A. 10) yields the complete cubic variation, 



ab 



(A.15) 



StotH = 



J d 3 x \{-di{q ln (d a q nb + d b q na - d n q ab )) - \d a q m id b q ml - d m q a id l q bm ^P ab 



+ \tr{qqP) + \q r b n Q ab p a 



(A.16) 



In order to cancel this variation, as already explained in text, we may only modify 5^p by adding 
to it terms proportional to q 2 , while preserving the symplectic form without help from other sectors. 
It produces, in the variation of the symplectic form, a contribution of the type q 2 q that cannot 
originate from anywhere else. 

In this remaining variation (A.16), we see two "non-interacting" sectors: the first is non-local 
(V~ 2 ) and the second local; this splitting is unambiguous. Inspecting the position of the indices of 
the derivatives in the non-local part, one sees that they cannot produce further local contribution 
by integration by parts. Hence the two contribution must be cancelled by separate terms in the 
variation S^p ab . Let us dispose first of the local sector. Here, we modify 5^p ab as follows 



5 m p ab = 0{qq)ab _ qacQcb + (6 <_> a) . 

The resulting variation of the symplectic form vanishes. 

S q 2S = J d A x (\q ab e mra d r (q mc q cb ) - \e mra d r (q mc q cb )q ab ) = 0. 



(A.17) 



(A.18) 



As explained in text, any term ~ J (Oq)qq effectively changes sign upon simultaneous space and 
time integrations by part. This brings us down to the non-local part, the final (and as we now see, 
incurable) obstacles: 



SnlHc = - J d 3 x di(q ln (d a q nb + d b q na - d n q ab )) + \ d a q„ d d b q ml + (d m q al d l q bm ) 



ab 



(A.19) 



This proliferation of terms represents the additional index possibilities once the extra V 2 brings 
in two extra derivatives in the numerator. 



Write (A.19) schematically as 

SnlH c = J d 3 x[qq]P = J d 3 xp6[pp} . 



(A.20) 



As in the local term's analysis, (A.20) can be cancelled by S NL p ab = {0[qq]} ab acting on the p 

(2) 

part of Hq, and the last hurdle is the ensuing symplectic variation, engendered by 5 NL p: 



J qqOq , 



(A.21) 
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which must again vanish by itself, there being no further variations left to help it. From (A. 19), 
then, 

6 NL S = J d± X {d em ( qa£qbm - Qabqem ) + \ qmMb + 2d e ( qemqrna , b ))d qab . (A.22) 

The terms in curly brackets are essentially all the possible forms involving qq and dd as a 2- 
index tensor: two, one and no "dummy" derivatives. Unlike in the local term's qqq there are no 
cancellations here and indeed there are concrete counterexamples of TT ^-tensors with SnlS / 0. 
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